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The blackbody radiation (BBR) shift of the 5s — 4^5/2 clock transition in ^^Sr""" is calculated to be 
0.250(9) Hz at room temperature, T = 300 K, using the relativistic all-order method where all single 
and double excitations of the Dirac-Fock wave function are included to all orders of perturbation 
theory. The BBR shift is a major component in the uncertainty budget of the optical frequency 
standard based on the ^^Sr"'" trapped ion. The scalar polarizabilities of the 5s and 4ci5/2 levels, as 
well as the tensor polarizability of the 4^5/2 level, are presented together with the evaluation of their 
uncertainties. The lifetimes of the 4ci3/2, 4d5/2, 5pi/2, and 5p3/2 states are calculated and compared 
with experimental values. 
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I. INTRODUCTION 

The current definition of a second in the International 
System of Units (SI) is based on the microwave transition 
between the two hyperfine levels of the ground state of 
^^•^Cs. The present relative standard uncertainty of Cs 
microwave frequency standard is around 4 x 10~^^ 
Since the frequencies of feasible optical clock transitions 
are five orders of magnitude greater than the standard 
microwave transitions, optical frequency standards may 
achieve even smaller relative uncertainties. Significant 
recent progress in optical spectroscopy and measurement 
techniques has led to the achievement of relative stan- 
dard uncertainties in optical frequency standards that 
are comparable to the Cs microwave benchmark. 

In 2006, the International Committee for Weights and 
Measures (CIPM) recommended that the following tran- 
sitions frequencies shall be used as secondary representa- 
tions of the second i2ll : ground-state hyperfine microwave 
3], 5s ^5*1/2 — 4d ^-D5/2 optical 



transition in ®^Rb ^. 

transition of the ^^St+ ion 0, % h(f^Qs ^Si/2{F = 
0) - bdP&s'^ '^D5/2{F = 2) optical transition in ^^^Rg+ 
ion 0, i, 6s ^Si/2{F = 0) - 5d "^D^/^iF = 2) optical 
transition in ^''^Yh+ ion d, ^ and Ss^ \SVi - 5s5p '^Pq 
transition in '^'^Sr neutral atom [H [H, [l3|. With ex- 
tremely low systematic perturbations and better stabil- 
ity and accuracy, such optical frequency standards can 
reach a systematic fractional uncertainty of the order of 
10-18 [|,|l|. More precise frequency standards will open 
ways to improve global positioning systems and tracking 
of deep-space probes, and perform more accurate mea- 
surements of the fundamental constants and testing of 
physics postulates. 
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In this paper, we treat one of the optical transitions 
recommended as secondary standard: the 5s ^6*1/2 — 
Ad electric-quadrupole transition of ^^Sr"^ at 445 

THz (674 nm). The reported frequency measurements of 
this transition in a single trapped ^^Sr"*" ion have achieved 
spectral resolution of better than 1.5 Hz 0, 0, [l^, [l^ . 
Methods based on this transition have the potential to re- 
duce relative systematic uncertainty to the level of 10" 
or below [6|. 

The accuracy of optical frequency standards is limited 
by the frequency shift in the clock transition caused by 
the interactions of the ion with external fields. The lin- 
ear Zeeman shift in the ^^Sr"*" system can be eliminated 
by use of the line center of symmetric Zeeman states; the 
second-order Zeeman shift is around 1 mHz and is neg- 
ligible at current level of precision Q . The second order 
Doppler shifts due to micromotion of the trapped ion are 
estimated to be less than 0.01 Hz [61]. The major con- 
tributions to the systematic frequency shifts come from 
Stark shifts with the blackbody radiation (BBR) Stark 
shift being one of the most important contributions at 
room temperature. Experimental measurements of the 
BBR radiation shifts are difficult. Here, we present the- 
oretical calculations that result in estimates of this shift 
that reduce the previous uncertainty [S, by a factor of 
10. 

In this paper, we present a relativistic all-order calcu- 
lation of the static polarizabilities of the 5s 1/2 and 4^5/2 
states of *8Sr+. The relativistic all-order method used 
here is one of the most accurate methods used for the 
calculation of atomic properties of monovalent systems. 
Readers are referred to Ref. flT'l and references therein 
for a review of this method and its applications. We 
use these polarizability values to evaluate the BBR shift 
of the clock transition at room temperature. The dy- 
namic correction to the electric-dipole contribution and 
multipolar corrections due to Ml and E2 transitions are 
incorporated. The uncertainty of the final BBR shift is 
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estimated to be 3.6%. Lifetimes of the low-lying excited 
4^3/2, 4(i5/2, 5pi/2, and 5p3/2 states are also calculated 
and compared with experiments. 



II. METHOD 

The electrical field E radiated by a blackbody at tem- 
perature T, as given by Planck's law, 



E'^(uj)duj 



8a3 



LO duj 



TT exp(a;/fcBr) — 1 ' 



(1) 



induces a nonresonant perturbation of the optical transi- 
tion at room temperature . Assuming that the system 
evolves adiabatically, the frequency shift of an atomic 
state due to such an electrical field can be related to the 
static electric-dipole polarizability by (see Ref. [l9j |) 
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(831.9 V/m)- 



(T{K)y 
V 300 / 



ao(l + r;). (2) 



Here, 77 is a small dynamic correction due to the fre- 
quency distribution. Only the electric-dipole transition 
part of the contribution is considered in the formula be- 
cause the contributions from Ml and E2 transitions are 
suppressed by a factor of [l^ . We estimate these mul- 
tipolar corrections together with the dynamic correction 
ri in Section |V] of this work. The overall BBR shift of 
the clock transition frequency is then calculated as the 
difference between the BBR shifts of the individual levels 
involved in the transition: 

Abbr(5s ^ 4^5/2) = -^[ao(4d5/2) - ao(5si/2)] 

X (831.9V/m)2(^ffi)y. (3) 

Therefore, the evaluation of the BBR shift requires ac- 
curate calculation of static scalar polarizabilities of the 
®*Sr+ in the 5si/2 ground and 4^5/2 excited states. The 
static tensor polarizability of the 4^5/2 state is also calcu- 
lated in the present work though the effect of the tensor 
part of polarizability is averaged out due to the isotropic 
nature of the electric field radiated by the blackbody. 

The calculation of the scalar polarizability of a mono- 
valent atom can be separated into three parts: the con- 
tribution of the electrons in the ionic core, cicorc; a small 
term, a^c, that changes the core polarizability due to the 
presence of the valence electron; and the dominant con- 
tribution, a^, from the valence electron. The ionic core 
polarizability used here was calculated using the random- 
phase approximation (RPA) We calculate the a^c 
contribution in the RPA approximation as well for con- 
sistency with the ionic core value. The valence scalar uq 
and tensor a2 polarizabilities of an atom in a state v can 
be expressed as the sum over all intermediate states k 



TABLE I: Contributions to the 5s 51/2 scalar (qq) static 
polarizabilities in ^^Sr"*" and their uncertainties in units of Aq- 
The values of corresponding El matrix elements are given in 
eao. 



Contribution 


(fcPI|5si/2> 


do 


5Si/2 - 5pi/2 


3.078 


29.23 


5Si/2 - 6pi/2 


0.025 


0.001 


OSi/2 - 7pi/2 


U.OOO 


0.004 


5Si/2 - 8pi/2 


0.054 


0.003 


5Si/2 - 5p3/2 


4.351 


56.48 


5Si/2 - 6p3/2 


0.034 


0.002 


5Si/2 - 7p3/2 


0.053 


0.003 


5Si/2 - 8p3/2 


0.054 


0.003 


dcorc 




5.81 






-0.26 


Otail 




0.02 


Ototal 




91.30 



allowed by the electric-dipole selection rules: 



3(2j„ 



Ek — Ey 



(4) 



a2 = -4C^(-1) 



ifjv 1 Jk \ {k\\D\\v) 



C = 



5j.(2j„ - 1) 



1 > 2 / Ek-E, 
1/2 



6(j,-f 2)(2j;, + l)(2j,„ + 3) 



(5) 



where (fc||D||u) are the reduced electric-dipole (El) ma- 
trix elements and Ei is the energy of the ith state. We 
also separate the valence polarizability into two parts, 
the main term amain containing the first few dominant 
contributions and the remainder ataii- We use electric- 
dipole matrix elements calculated using the relativistic 
single-double (SD) all-order method (se e jlTj for detail) 
and experimental energies from Ref. [2l| for the calcula- 
tions of the main term. Triple excitations are included 
partially where needed; the resulting values are referred 
to as SDpT (single, double, partial triple) data. 

The tail contribution for the 5s state is negligible and is 
estimated in the lowest-order Dirac-Fock (DF) approxi- 
mation. Significantly larger tail contribution to the 4(i5/2 
polarizability is evaluated in both Dirac-Fock (DF) and 
RPA approximation and scaled to account for other miss- 
ing correlation corrections. 

In this work, we use atomic units (a.u.), in which, e, 
TOe, 47reo and the reduced Planck constant h have the nu- 
merical value 1. Polarizability in a.u. has the dimension 
of volume, and its numerical values presented here are 
thus expressed in units of ajj, where ag w 0.052918 nm is 
the Bohr radius. The atomic units for a can be converted 
to SI units via a / h[}iz / {V / m^] = 2.48832 x 10'^a[a.u.], 
where the conversion coefficient is AneoaQ/h and Planck 
constant h is factored out. 

We have used the B-spline method to construct a fi- 
nite basis set for radial Dirac equations as introduced in 
Ref. [l^l. 70 B-splines of order k = 8 are constrained to 
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TABLE II: Contributions to the 4dr,/2 scalar (ao) and tensor 
(a2) static polarizabilities in ^'^Sr'^ and their uncertainties in 
units of Oq. The absolute values of corresponding El reduced 
matrix elements are given in eao. 



Contribution 


\K II 1^11^(25/2/ 




(y.2 


A ^ 

4a5/2 


- 5p3/2 


■4:. 10 ( 


^4:. ±U(^zy J 


A_A 1 f<(9Q\ 

-4:4:. LUy^iJ J 


A ^ 

4a5/2 


- 6p3/2 


n 1 /19 

U. l^Z 


01 9/'9'l 


01 9^'9^ 


A ^ 

4a5/2 


'P3/2 


U.U / 


00"^ 


00"^ 

-u.uuo 


A ^ 

405/2 


- OP3/2 


U.UtJO 


001 


001 


A ^ 

4a5/2 


/I f 

- 4/5/2 






u.o ( uy^ J 


Afl 


- 0/5/2 


0.442 


os^/9'1 


Q7(9^ 
u. y 1 \^ J 


Afl 

405/2 


- 0,15/2 


n 9Q7 




040 
u.u^u 


A ^ 

4055/2 


- '75/2 


910 


018 


091 

U.UZ 1 


A ^ 

4o!5/2 


- 0/5/2 


1 "^7 


OOQ 


01 

U.U -LU 


A^ 

40!5/2 


- y/5/2 


1 "^S 
U. loo 


007 
U.UU i 


008 
U.UUo 


Afl 

40!5/2 


1 n f 

— J^U/5/2 




00^ 


00^0 

U.UUiJU 


A ^ 

40!5/2 


- 11/5/2 


OQS 


00"^ 

U.UUO 


004 
u.uu^ 


Art , 


— t^/5/2 


ORf^ 


00*^ 

U.UUO 


00"^ 

U.UUO 


A^ 

40!5/2 


A f 
- 4/7/2 




fi ^7fi/'70'l 


-Zr. 0^4:01^^)0 J 


4^5/2 


— 5/7/2 


1 Q7Q 


1 .uyy i^ou J 


f^07r 1 1 
-U.UU t L J 


4d5/2 


- 6/7/2 


1.329 


0.698(11) 


-0.249(4) 


4d5/2 


— 7/7/2 


0.979 


0.360(5) 


-0.128(2) 


4d5/2 


— 8/7/2 


0.764 


0.212(4) 


-0.076(2) 


4d5/2 


- 9/7/2 


0.619 


0.136(2) 


-0.049(1) 


4d5/2 


- 10/7/2 


0.517 


0.093(1) 


-0.033 


4d5/2 


- 11/7/2 


0.440 


0.067(1) 


-0.024 


4d5/2 


- 12/7/2 


0.381 


0.050(1) 


-0.018 








5.81(29) 




Ovc 






-0.40(10) 




Otail 






2.06(20) 


-0.59(7) 


Qtotal 






62.0(5) 


-47.7(3) 



a spherical cavity with R = 220 a.u. for each angular 
momentum. Such a large cavity is chosen to accurately 
evaluate as many 4^5 /2 — n/7 /2 transitions as practically 
possible to reduce the uncertainty in the remainder. 



III. POLARIZABILITIES 

Table [J shows the contributions of the individual tran- 
sitions to the ground state scalar polarizability ao- The 
main contributions are listed separately along with the 
respective values of the electric-dipole matrix elements. 
The tail contributions are grouped together as ataii- For 
the main contributions, we use our ab initio SD all-order 
values of the matrix elements and experimental energies 
from Ref. [2l[. The 5si/2 — 5pi/2 and 5si/2 — 5p3/2 tran- 
sitions contribute over 99.9% to the valence polarizabil- 
ity and 94% of the total polarizability value. The same 
calculation of these transitions in Rb agrees with high- 
precision experiment to 0.26% [2^. In fact, the SD val- 
ues for the primary ns — np transitions in Li, Na, K, Rb, 
and Cs agree with various types of high-precision exper- 
iments to 0.1% - 0.4% [2I]. There is no reason to expect 
reduced accuracy in the case of Sr"*", and we take the 
uncertainty of these matrix element values to be 0.5%. 
Unfortunately, we know of no way to accurately estimate 



the missing additional contributions to the dominant cor- 
relation correction to these transitions, unlike the case of 
the Ad — bp and 5s — Ad transitions, where semiempirical 
scaling makes possible an uncertainty estimate that does 
not directly depend upon comparison with experiment. 

The core contribution taken from Ref. [l^l is estimated 
to be accurate to 5%, based on the comparison of the 
RPA and experimental polarizability values for noble 
gases. The tail contribution is calculated using the DF 
approximation and is negligible in comparison with the 
total polarizability. The error of the tail is taken to be 
100%. As a result, all uncertainties except the ones as- 
sociated with the 5s — 5p matrix elements are negligible. 
The resulting final uncertainty of the 5s polarizability is 
thus estimated to be 1%. We note that accurate measure- 
ment of either 5s — 5p oscillator strengths or bp lifetimes 
{bp — M contributions are small and can be accurately 
calculated) will help to significantly reduce this uncer- 
tainty. 

Table shows the contributions from the individual 
transitions to the 4^5/2 polarizability. Three types of 
transitions contribute to the 4c?5/2 polarizability: 4^5/2 — 
np3 /2 , 4(^5 /2 - n/5/2 , and id^ /2 - nf^ /2 . The sum over the 
4*^5/2 ~ transitions converges very quickly with the 

4^5/2 ~ 5p3/2 term being overwhelmingly dominant. We 
obtain an accurate value for this matrix element using 
a semi-empirical scaling procedure that evaluates some 
classes of correlation corrections omitted by the current 
all-order calculations. The scaling procedure is described 
in Refs. [23,[2^. Briefly, the single valence excitation co- 
efficients are multiplied by the ratio of the correspond- 
ing experimental and theoretical correlation energies, and 
the matrix element calculation is repeated with the mod- 
ified excitation coefficients. The scaling procedure is par- 
ticularly suitable for this transition because the matrix 
element contribution containing the single valence excita- 
tion coefficients is dominant in this case (but not for the 
5s — bp matrix elements discussed earlier). We conduct 
the scaling starting from both SD and SDpT approxima- 
tions. The scaling factors for the SD and SDpT calcu- 
lations are different, and we take scaled SD value as the 
final result for the 4^5/2 — 5p3/2 matrix element, based 
on the comparisons of similar calculations in alkali-metal 
atoms with experiments [H, [13, [H, HI- The absolute 
values of the reduced 4^5/2 — 5p3/2 matrix elements cal- 
culated in different approximations are summarized in 
Table Hill together with four other transitions that rep- 
resent similar cases. The uncertainties are determined 
as the maximum difference between the scaled SD values 
and the ab initio SDpT and scaled SDpT values. A no- 
table feature of this table is close agreement of the scaled 
SD and SDpT results. The sum of the contributions from 
the 4(i5/2 —nf^/2 and 4^5/2 — 're/7/2 transitions converges 
slowly; therefore we include as many transitions as re- 
alistically possible in the main term calculation. Scaled 
values are used for the 4^5/2 — 4/5/2 and 4^5/2 — 4/7/2 
transitions in the polarizability calculations. 

The tail contribution of the 4c?5 /2 — n/7 /2 terms is par- 
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TABLE III: Reduced electric-dipole transition matrix elements calculated using different approximations: Dirac-Fock (DF), 
single-double all-order method (SD), and single-double all-order method including partial triple-excitation contributions 
(SDpT); the label "sc" indicates the corresponding scaled values. All values are given in atomic units. 



Transition 


DF 


SD 


SDpT 


SDec 


SDpTsc 


Final 


4ds/2 - 5p3/2 


5.002 


4.150 


4.198 


4.187 


4.173 


4.187(14) 


4^5/2 - 4/5/2 


0.964 


0.779 


0.790 


0.789 


0.785 


0.789(4) 


4^5/2 - 4/7/2 


4.313 


3.486 


3.536 


3.528 


3.509 


3.528(19) 


4^3/2 - 5pi/2 


3.729 


3.083 


3.119 


3.112 


3.102 


3.112(10) 


4^3/2 - 5p3/2 


1.657 


1.369 


1.386 


1.383 


1.378 


1.383(5) 



TABLE IV: Contributions to the lifetimes of the 5pi/2 and 
5p3/2 states. The transitions rates A are given in 10"^ s"^ and 
the lifetimes r are given in ns. 



TABLE V: Lifetimes of the 4^3/2 and 4^5/2 levels (s) 



5pi/2 



A(5pi/2 - 5s) 128.04 

A(5pi/2 - 4^3/2) 7.54 
135.58 

t(5Pi/2) 

Present 7.376 

Expt. [30] 7.39(7) 

Expt. [31] 7.47(7) 



5p3/2 



^(5p3/2 - 5s) 141.29 

A(5p3/2- 4^3/2) 0.96 

A(5p3/2 -4d5/2) 8.06 

J2A 150.31 

t(5P3/2) 

Present 6.653 

Expt. [30] 6.63(7) 

Expt. [31] 6.69(7) 



ticularly large; its DF value (3.5 a.u.) is 5% of the total 
polarizability. Therefore, we carry out several additional 
calculations to accurately evaluate the tail contribution 
and estimate its uncertainty. Since the largest part of 
the correlation correction for the 4c?5/2 — ?^/7/2 transi- 
tions with n > 9 terms comes from RPA-like terms, the 
RPA approximation is expected to produce a better re- 
sult than the DF one. We carried out the RPA calcu- 
lation of the tail and obtained a lower value of 2.9 a.u. 
We also calculated the last few main terms using the DF 
and RPA approximations and compared the results with 
our all-order values. We found that the DF and RPA 
approximations overestimate the polarizability contribu- 
tions by 35% and 28%, relative to DF and RPA values. 
To improve our accuracy, we scale both DF and RPA re- 
sults by these respective amounts to obtain a DF-scaled 
value of 2.26 a.u. and RPA-scaled value of 2.06 a.u. We 
take RPA-scaled value as the final one and the difference 
of these two values as its uncertainty. 

We also list the contributions from various transitions 
to the 4(^5 /2 tensor polarizability 02 in Table [ill The 
4^5/2 — np2/2 transition gives the dominant contribu- 
tion to the tensor polarizability. The tail contribution 
is smaller yet significant and is obtained by the same 
procedure as the tail of the scalar 4^5/2 polarizability. 



IV. LIFETIMES 



Levels Experiment 



Other calculations Present 



4d3/2 0.435(4) [32] 


0.443 [33] 


0.441(3) 


0.455(29) \32\ 


0.426(8) \3£ 




0.435(4) [35] 


0.422 [32] 






0.441 [36] 




4d5/2 0.3908(16) [37] 


0.404 [33] 


0.394(3) 


0.408(22) [32] 


0.357(12) [34] 




0.372(25) [38] 


0.384 [32] 






0.396 [36] 





The contributions to the lifetimes of the 5pi/2 and 
5p3/2 states are given in Table HVl Experimental ener- 



gies from Ref. [2l|] are used in the evaluation of the tran- 
sition rates. The lifetime is calculated as the inverse of 
the sum of the appropriate Einstein A-coefficients, which 
are proportional to the square of the dipole matrix el- 
ements; experimental energies from Ref. [2l| are used. 
Our results are in excellent agreement with experimental 
lifetimes r(5pi/2) — 7.39(7) ns and T(5p3/2) — 6.63(7) 
ns by Pinnington et al. measured using a laser-induced 
fluorescence method [13, and t(5pi/2) = 7.47(7) ns and 
t(5p3/2) = 6.69(7) ns by Kuske et.al. measured using 
fast-beam-laser technique ^31i] . 

As a further test of accuracy of our approach and 
accuracy of our all-order Ad wave functions, we carry 
out the calculation of the 4c?3/2 and 4^5/2 lifetimes 
that requires evaluation of the electric-quadrupole and 
magnetic-dipole transitions. 

The lifetime of the 4c?3/2 state is calculated to be 
0.441(3) s, where the main contribution comes from the 
i?2(4(i3/2 — 5S1/2) = 11.13(3) a.u. matrix element. The 
contribution from 4^3/2— 5si/2 Ml transition is evaluated 
to be negligible. The most recent lifetime measurements 
of 4(i3/2 states of Sr+ include 0.435(4) s result obtained 
by using optical pumping and 0.455(29) s value obtained 
by using laser probing as reported in the same work [s^ . 
Our result agrees with the experimental values within 
the uncertainty limits. The lifetime of the 4^5/2 state 
is calculated to be 0.394(3) s. The contribution to the 
A-coefficients from the 4^5/2 — 5s 1/2 E2 transition over- 
whelmingly dominates, and the corresponding reduced 
matrix element is 13.75(4) a.u. The reduced matrix ele- 
ment for the 4:d^/2 — 4^3/2 E2 transition is 5.98(2) a.u., 
but its contribution to the lifetime is negligible due to 
the small energy interval between these two states. Our 
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TABLE VI: Comparison of static scalar polarizabilities for the 5si/2 and 4^5/2 states and blackbody radiation shift for the 
5si/2 — 4d5/2 transition in *^*Sr^ ion at T = 300 K. The polarizabihty values are in ajj and BBR shift is in Hz. 





Present work 


Ref. [33] 


Ref. [39] 


Ref. [5] 


Ref. [40] 


Ref [6] 


Qo(5si/2) 


91.3(9) 


89.88 


93.3 


84.6(3.6) 


91.47 




00(4^5/2) 


62.0(5) 


61.77 62.92* 


57.0 


48(12) 






BBR shift 


0.250(9) 


0.242 0.233* 


0.31 


0.33(0.12) 




0.33(9) 



*Rcsults are obtained using experimental energies. 



value for the 4c?5/2 — 4c?3/2 Ml reduced matrix element, 
1.55 a.u., is in agreement with the result from Ref. [3^ . 
The contributions from Ml transitions only affect the 5th 
decimal of the lifetime result of the 4^5/2 state, and can 
be neglected at the present level of accuracy. The two- 
photon transitions contribute 0.03% to the Ad lifetimes. 
Our result is found in good agreement with the lifetime 
measurements of 0.3908(16) s and 0.408(22) s done by 
Letchumanan et al. [s^ and Biemont et al. [13], respec- 
tively. More experimental and theoretical results for the 
lifetimes of these two states are given in Table IVl 

V. BBR SHIFT 

We use our scalar polarizability values to evaluate the 
shift in the clock transition in ^^Sr+ due to blackbody 
radiation at T = 300 K to be 0.252(9) Hz. The dy- 
namic correction [l^ is estimated to be ?7=0.0013 and 
7^=0.0064 for the 5s and 4^3/2 states, respectively. The 
resulting correction to the BBR shift is —0.002 Hz and 
our final value for the BBR shift is 0.250(9) Hz. The 
overall uncertainty in the final result comes from the un- 
certainty in the values of the 5s — 5p3/2 and — 4^5/2 
matrix elements, and 4^5/2 — tail- The first two 

sources of the uncertainties may be removed if these val- 
ues were determined experimentally. The 5s — 5p3/2 and 
5p3/2 ~ 4(^5/2 matrix elements can be obtained via either 
lifetime, ground state polarizability, oscillator strength, 
or light shift ratio measurements [4l|, with the first two 
type of experiments useful for the first matrix element 
and second two types for both matrix elements. We note 
that ionic core uncertainty is not included in the BBR 
uncertainty since the core contribution is the same for 
both levels and subtracts out. We note that the small 
term a^c that changes the core polarizability due to the 
presence of the valence electron is different for the 5s and 
4(^5/2 states and contributes 0.5% to the BBR shift. We 
also estimated the Ml and E2 contributions to the BBR 
shift using the approach described in Ref. [I9j and found 
them to be negligible (below 0.01%). 

In Table IVIi we compare our polarizability and BBR 
shift results with other theoretical calculations [1, 0, 
m, [13, [13]. We note that our calculation is the most 
complete one at present time. The calculation of Mitroy 
et al. [33! carried out by diagonalizing a semi-empirical 



Hamiltonian in a large-dimension single-electron basis 
gives the BBR shift of 0.242 Hz. The use of experimen- 
tal energies changes this result to 0.233 Hz [Sj]- The 
accuracy of the scalar polarizability of the ground state 
ao = 89.88 was estimated to be 2 - 3% in Ref. [11]. The 
results of Barklem and O'Mara [13] are derived from ex- 
perimental oscillator strengths (note that the core po- 
larizability of 5.8ag was added to the values listed in 
Ref. [13). The results of Madej et.al. [| are obtained 
mainly by summing over the transition rates calculated 
in Ref. [42| using the multiconfiguration Hartree-Fock 
(MCHF) method. It appears that the core contribution 
(5.8 a.u.) was omitted in [3]. Addition of the core contri- 
bution to the results of [1] leads to ao(5s) = 90.4 a.u. and 
00(4^5/2) = 54 a.u. values. Omission of the higher-order 
transition contribution resulted in lower values for the 
scalar polarizability of the 4(i5/2 state in [^. The scalar 
polarizability of the 5s 1/2 state calculated by Patil and 
Tang [40] is obtained by evaluating the transition ma- 
trix elements with simple wave functions based on the 
asymptotic behavior and on the binding energies of the 
valence electron. Our result is in good agreement with 
their calculation. Another estimation of the BBR shift is 
given by Margolis et al. in Ref. (gI], but the approach is 
not stated. 



VI. CONCLUSION 

In summary, we calculated the polarizabilities of the 
5s and 4^5/2 states in ^^Sr+ and the value of BBR shift 
of the corresponding clock transition at room tempera- 
ture. The dynamic correction to the electric-dipole con- 
tribution and the multipolar corrections due to Ml and 
E2 transitions were estimated and found to be small at 
the present level of accuracy. Lifetimes of the low-lying 
excited 4^3/2, 4ci5/2, 5pi/2, and 5p3/2 states were also 
calculated and compared with experiments for further 
tests of our approach. The uncertainty of the final BBR 
value was estimated. The main contributions to the un- 
certainties were analyzed and possible experiments were 
suggested to further reduce the uncertainties of the BBR 
shift. 

This research was performed under the sponsorship of 
the US Department of Commerce, National Institute of 
Standards and Technology. 
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